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ABSTRACT 

The standard Hartree-Fock approximation of the O(N)— invariant 4 model suffers 
from serious renormalization problems. In addition, when the symmetry is spontaneously 
broken, another shortcoming appears in relation to the Goldstone bosons: they fail to 
be massless in the intermediate states. In this work, within the framework of out-of- 
equilibrium Quantum Field Theory, we propose a class of systematic improvements of the 
Hartree-Fock resummation which overcomes all the above mentioned difficulties while 
ensuring also exact Renormalization-Group invariance to one loop. 

PACS: ll.10.Ef, ll.10.Gh, ll.10.Pg 

1 Introduction, summary and outlook 

The research in out-of-equilibrium dynamics of quantum fields has received, in recent 
years, a great impulse by cosmology as well as by particle and condensed matter physics. 
Indeed, a first-principle theoretical treatment of many important phenomena, such as 
the reheating of the universe after inflation or the thermalization of the quark gluon 
plasma in the ultra-relativistic heavy-ion colliders (RHIC, LHC), is required for a good 
qualitative and quantitative understanding of the late time and strongly coupled evolution 
of quantum field systems. 

In particular, this necessity has encouraged the study of nonperturbative approaches 
to Quantum Field Theory (QFT) that could provide all-orders partial resummations of 
Feynman diagrams PjJ |2] . In fact standard perturbation theory does not yield satisfactory 
results, except for very short times, when non-equilibrium conditions are involved. 

Mean-field approximations such as Leading-order large-N expansion [3J EJ E] and 
Hartree, or Hartree-Fock (HF) variational method |7J |HJ HH EI] are the simplest and 
most studied fTJ El El El El Ej resummation schemes. Their main features are well 
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known: they do provide a backreaction term on the evolution of quantum fluctuations 
that stabilize dynamics after parametric amplifications or spinodal instabilities, but they 
fail to reproduce an important property of late time dynamics such as thermalization. 

A more powerful tool is the 2PI (or 2PPI) effective action [T71 ITS] whose expansions at 
two (or more) loops order or at next-to-leading order in 1/N provide resummations that 
go beyond mean field approaches [2111 El 1221 ES] and yield, indeed, approximate numerical 
thermalization at strong coupling. 

More formal aspects of resummed approximations, such as their renormalization prop- 
erties, have been studied as well |2H1 |23j . Recently a systematic method has been put 
forward [23 EDI 123 EUj that removes divergences in the ^-derivable approximations by 
applying a BPHZ subtraction procedure to diagrams with resummed propagators. 

In j20j we considered the simple HF approximation of the O(N) 4 model in the 
unbroken symmetry phase. It is known that it cannot be consistently renormalized by 
the usual renormalization of bare coupling and mass [SI [23 121] • We showed that this non- 
renormalizability is due to the absence of leading logarithmically divergent contributions 
coming from diagrams which do not have a "daisy" topology and therefore are not present 
in the standard HF resummation. The inclusion of these contribution together with 
suitably chosen finite parts led us to the definition of a renormalized and Renormalizat ion- 
Group invariant version of HF equations. 

In the present paper we further develop this analysis, making it more systematic and 
transparent, and then apply it to the subtler case of spontaneously broken O(N) sym- 
metry. In this case, together with the shortcomings as far as pure renormalizability is 
concerned, the HF approximation shows an unphysical non-gapless behaviour. Namely, 
the effective resummed propagators of the Goldstone bosons, as defined by the HF equa- 
tion of motion, are not massless. This is a well known problem which has been cured 
by defining suitably modified approximations [211 E2] • In this paper we define a system- 
atic improvement of the standard HF approach which ensures gaplessness together with 
renormalizability and RG-invariance. In practice we include all leading logarithmically 
divergent contributions needed to render finite the resummed propagator masses (which 
are log A dependent in standard HF) while choosing the finite parts in such a way to have 
massless transverse degrees of freedom. 

More in detail, we proceed as follows. First, by analysing the diagrammatic resumma- 
tion performed by the HF approximation, we identify the missing leading logarithmically 
divergent contributions that cause the non-renormalizability. As will be apparent below, 
this analysis closely parallels that in [30J since, after all, the UV behavior of the the- 
ory is the same regardless of the breaking of the O(N) symmetry. Then we construct 
our modified HF approximation, including the missing contributions, while choosing the 
corresponding finite parts by means of the following recipe 

- determine the general features of the HF approach that must be shared by the modi- 
fied approximation, such as mean-field structure, O(N) Ward identities, leading-log 
structure and more, and parameterize the class of resummations having such prop- 
erties. 

- identify amid this family of approximations the one(s) with suitable properties of 



2 



renormalizability, RG-invariance, gaplessness and infrared finiteness. 

One important difference from [30J is that, when the O(N) symmetry is spontaneously 
broken, the above procedure does not lead to a unique result except for the N = 1 case 
(Z 2 symmetry). When N > 1 we find a whole class of approximations, corresponding to 
different choices of finite parts, that share the required characteristics. Like in (3U| the 
improved equations have, by construction, a mean field structure but, unlike the standard 
ones, are nonlocal in space and time. 

In sec. 121 we review some general concepts. We introduce the CTP formulation of 
out-of-equilibrium problems and define the HF approximation as resummation of bubble 
diagrams recalling the corresponding, well known, equations of motion [see eqs. (|2.4|) and 
eqs. (12. 6|) ]. Some important general features of the mean field approximations are pointed 
out. 

In sec. El we study, as an illuminating example, the case of spontaneously broken Z 2 
symmetry. Subsec. 13.11 is dedicated to a general analysis of the approximation. In sub- 
sec. HOI we apply the standard renormalization procedure and point out its shortcomings. 
In subsec. 13. 31 we define our modified HF approximation which is renormalizable and RG- 
invariant. As we already said in this case we are led to a unique result [see eqs. fl3.38|) 
and eqs. 1)3.42)1 ] . Further remarks on the freedom of choosing various initial conditions 
without spoiling renormalization and RG-invariance are made in subsec. 13.41 

In sec. |l]we generalize to N > 1, that is to a continuous O(N) spontaneous symmetry 
breaking, by first analysing some general features [subsec. l4.1j of the standard HF approx- 
imation and its renormalization properties [subsec. 14.2) and then defining our modified 
approximation by the general recipe given above [subsec. 14.3) . The final results will be a 
set of constraints that should be satisfied by renormalizable and RG-invariant Hartree- 
Fock-like approximations together with two examples [see eqs. (|4.38|) and eqs. (|4.39|) ]. 

There are several possible developments along the lines of this work. First of all it 
would be interesting to study other properties of the different improved approximations 
that we have constructed in the N > 1 case and, hopefully, identify further constraints 
that would allow to single out a unique result. Secondly, a numerical study of the modified 
HF field equations found in this paper and in [3U] would be needed to investigate how 
the space-time nonlocalities affect the time evolution as compared to the standard HF 
approximation, which is known to fail even qualitatively at late times. Another challenging 
task is the extension of our approach to the full two-loop 2PI effective action, through the 
inclusion of the nonlocal sunset diagram which is absent by definition in any mean-field 
approximation. In fact one should expect that, even if such an inclusion allows to recover 
renormalizability as compared to the conventional HF approximation, the two-loop 2PI 
self-consistent equations still lack RG invariance with the two-loop beta function, since 
the 2PI effective action does not contain all the diagrams which contribute to the next- 
to-leading ultraviolet divergences. 
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2 Generalities 



To study the real-time dynamics from a given set of initial conditions we need to evaluate 
matrix elements of the form 

(*(t0)||O(tl,t2,...)l*(*0)> 

where is a generic state prepared at a initial time to and O is some operator depend- 
ing on the field evaluated at times U > to- Out-of-equilibrium QFT provides the general 
setup for the calculation such matrix elements, as well as more general expectation val- 
ues in statistical mixtures of pure states such as |\&). In this section we briefly review 
some generalities about non-equilibrium QFT, define the Hartree-Fock approximation 
and point out some properties that will be useful to derive the forthcoming results. In 
doing this, we restrict to the case of interest for this paper: scalar field theory in 3 + 1 
dimensions with quartic interaction and spontaneously broken O(N) symmetry. 

The field variables are <fi(x) where x = (xi,X2,Xs) are space coordinates and i = 
1,...,N is the O(N) index. Classical dynamics is defined by an action, functional of 
trajectories <Pi(x), with x = (x,t), in the form 



f a f 1 m 2 A 

S[<p] = d 4 x< -d^i{x)d^^i{x) - —cpi(x)(pi(x) - —\ Vi (x) Vi (x)f 



where m? and A are the squared mass (negative in the case of spontaneously broken 
symmetry ) and the coupling constant respectively. 

In QFT, the general approach to non-equilibrium dynamics was developed by Keldish 
and Schwinger and is known as closed time path (CTP) formalism. It allows to use 
standard functional methods (see jT7l [THJ ESI ESI)- by introducing path integrals on a 
time path going from t — to t — +oo and back. Field integration variables are then 
doubled and subdivided into (+)— components, for the path integral forward in time, and 
(—)— components for the backward piece. Given an initial state defined by the functional 
\I/ of the field configurations f(x), one writes down the functional integral 

e iw\j+,j-] = j T> (p+ 'D<p_itf[(p + ]y [tp_] e iS[<p+]-iS[<p-}+i{j + \<e+)-i(j-\<p-) p.!) 

where we have used the short-hand notation 

(a\M\b) = [ d 3 xd 3 y a j (x)M jk (x,y) b k (y) 



which will be useful later on. Integration in eq. (|2.1j) is on trajectories from t — to 
t = +oo (with the condition (p + = <p_ at t = +oo) and <p± in the wave functional is 
the t = section of (p±. Notice that, in the action functional S that enters in the path 
integral in eq. (j2.1|) . the constant parameters A, m 2 are to be substituted with the bare 
(cut-off dependent) ones A , rn%. In fact the theory should be thought as regularized with 
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an UV cut-off A and then renormalized to remove the divergent dependence on A. In the 
general setting, field renormalization should be included as well. Here we ignore it since, 
as we will see, it is absent in the HF approximation of the scalar theory. 

By construction W[j + , j_] is the generating functional of connected Green functions 



G_i — | (xi, . . . , x n , yi, . . . , y n ) 



-i 



,n+m 



Sj+(xi) . . . 5j+(x n )5j-.(yi) . . . 5j-(y m ) 
-i (tf 1 7{(p( yi ) . . . v{y m )}T {(p(xi) . . . (p(x n )} |*} c 



j+=o 

j-=0 



where T and T define time ordered and inverse ordered products, respectively, and in- 
ternal indices have been omitted for ease of notation. The effective action ripi[0 + ,0_], 
which is the generator of 1PI vertex functions, is the Legendre transform of W[j+,j-] 
from the currents j± to the fields <fi±. The equation of motion for the background field 
<f)(x) = (^| p(x) \ty) then reads 



5T 



ipi 



(2.2) 



Notice that the functional W as well as the IPI effective action parametrically depend on 
the initial state In our present discussion we consider an initial wave functional having 
the following Gaussian form 



*[<p] = Af exp [i(0(O)| tp - 0(0)) -(cp- 0(O)| [{g- 1 + iS] \cp - 0(0)) 

whose free parameters are the t = background field 0(cc,O), the t = background 
momentum <p(x, 0), the real symmetric positive kernel Qij(x,y) and the real symmetric 
kernel Sij(x, y). We can see that Tipi now depends parametrically only on the kernels Q 
and iS, while the t — background fields 0(0) and 0(0) enter instead as initial conditions 
for the equation eq. ()2.2|) . 

The perturbative diagrammatic expansion in the CTP formalism proceeds as in vac- 
uum QFT. The bare propagators are 



G { °i(x,y) = Gf(x,y) = -i (¥| T<p{x)<p{y) |¥) t 





[x, 


y) 


- 

= trp 


(b) 


[x, 


y) 




w 


[x, 


y) 


= G® 



free 



[x,y) = -i (^| T <p(x)<p(y) |^) conn | frcc 



G%l{x,y) = G%(y } x) = -i (*| <p(y)<p{x) |*) c 



I free 



Here the notation |f ree indicates that, in accordance to the tree level of the theory in the 
broken symmetry phase, the Heisenberg expectation values have to be calculated for a 
free scalar field theory with mass matrix mf- = -XoViVj, where Vi is the nonvanishing 
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vacuum expectation value. The bare vertices are 




* ^0 T ijkm v m 



' * ^0 T ijkm v m 



However, in order to obtain sensible results in a generic out-of-equilibrium contexts or 
even at equilibrium with nonzero temperature, it is known that calculations should go 
beyond plain perturbation and perform (partial) resummations to all orders in the cou- 
pling constant. A very successful resummation method is provided by the introduction 
of the 2PI effective action (see [T7t ITS]). It is defined as the double Legendre transform 
of the W generating functional with respect to the usual current one-point j± and to the 
two-points current K a p(x,y) coupled through the term 



d x Id y K af3 (x,y)(p a (x)tp (x) 



where a, (5 = ±. r 2 pi is a functional of the classical fields <p a and of the propagators G a p. 
It yields two equations of motions 



ST 



2PI 



q <5r 2 pi 



S(j) Q (x) * ' 5G a p(x,y) 



(2.3) 



Here the notation |* indicates that, by their physical meaning, the (±)-component fields 
and propagators have to satisfy, on the solutions of motion, the following relations 

<j)-{x) = 4>-{x) = 4>(x) 
G F (x, y) = G+-(y, x)6(x - y ) + G + _(x, y)9(y Q - x ) 
Gp(x, y) = G+_(x, y)9(x - y ) + G + ^(y, x)6(y - x ) 

Hence the system eq. ()2.3|) reduces to two coupled equations for and G + _ only. Moreover, 
any initial Gaussian state may be absorbed in the t = term for the j a and K a p currents, 
so that, by the double Legendre transform the initial Gaussian state disappears from the 
effective action, but fixes the initial conditions on <p and G + _. The role of 0(0) and 0(0) 
is immediate, while for the kernels we have 

G+-(x,y)\ Xo=yo= o = G(x,y) 

d 

—G+^(x,y)\ xo=yo=0 = 2i [GS\{x,y) + |5 3 (a? - y) 

Given T 2 pi at a certain perturbative loop order, if we solve the second equation in eqs. (|2.3J) 
for a generic and substitute the result G[<j>] into the first one we obtain the background 
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equations of motion corresponding to a resummed diagrammatic approximation of the 
1PI effective action Tipi. 

In the present scalar theory I^pi has the general form 

r 2P i [0, G] = S[<j>] + *-Tr [logG] + *-Tr [G^G] + T 2 [0, G] 

where S is the complete classical action of the double time path (i.e. S = S + — S_). 
Traces are taken over all indices i, a and x. Gq 1 is the second derivative of the action 
in a background, T 2 is the sum of all vacuum 2PI diagrams with G propagators and 
vertices defined by the classical action in a background. To two loops level the diagrams 
contributing to the r 2 are the "8" and "sunset" diagrams 



Now we can introduce the Hartree-Fock approximation as obtained considering only the 
first contribution (i.e. the "8" graph) to r 2 . 

r 2 = |A [G 2 F (x,x)-G%(x,x)] 

Notice that "8" is the only 2PI diagram made of "product" of loops corresponding to a 
mean field contribution to the mass. In the 1PI framework this corresponds to a resum- 
mation of all vacuum 1PI diagrams with daisy and superdaisy topologies of the form 



By explicitly using this form of T 2 in eqs. ()2.3|) , setting G(x, y) = | (x, y) + G- + (y, x)] 
and observing that the antisymmetric combination decouples, one obtains 

{[□ + m„ + |A o 0/ c (x)0 fc (x)l Sij + ^\ T ijkm G km (x,x)} 4>j{x) = 

(2.4) 

{(□ + m 2 ) + |A o r ijfcm [0 fc (x)0 m (x) + G km (x,x)] } G i:j (x,y) = 
where the tensor r is defined as 

Tijkm ^{.SijSkm SikSjm SimSjk) (^"^) 

It is convenient to introduce an equivalent formulation of the eqs. ()2.4|) in terms of mode 
functions Uh a (k is the wave vector and a the O(N) polarization), that will be the one 
used throughout the paper. For simplicity let us suppose that the initial (t = 0) kernels 
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are translationally invariant (while the background field and momentum may be point 
dependent). We can then write 



Gij{x,y) 



d 3 k 



Qij(k) e 



ikx 



d 3 k 



(27T 



<%(fc) e 



ikx 



,3 



(2vr) 3 ^ 

Next, we introduce the t = mode functions by 

u k , ai (x,0) = [g(k) 1 / 2 ] m e tk - x , u fc , of (aj,0)= [-f^(fc)- 1 + 2 ^(a;, 0) 



Then, one can easily verify that, by the identification 

d 3 k 



Gij(x,y) = Re J 



(27T) 



(v) 



eqs. f!2.4|) are equivalent to the following equations of motion 

f d 3 p } 

[n + m + l\ <f> k (x)<f>k(x)]8ij + |A r iifcTO J -j—^ u P)bm (x)u P)bn (x) j (f)j(x) 



J3 + mfybij + !AoTjj„,.,, 



m (x)0„(x) + 



d 3 p 
(2^)3 



Uk,aj(%) = 



(2.6) 



These are the HF equations of motions that will be used in the rest of the paper. 



We conclude this section by introducing a different representation of the CTP formal- 
ism, known as the physical representation (see j2E]) which allows us to derive some results 
useful below. We introduce the field redefinitions (omitting again the internal indices to 
simplify notation) 

0A = 0+~0- , c =-(0+ + 0_) 

and write the 1PI effective action as a functional of these new fields, I^pi = r[0A, <Pc] ■ By 
calculating vertex functions, one then finds 



8 n T 



and 



,(xi) . . . S(j) c (x n ) 



0A=O 



5<j) c (xi) . . . 5<fi c (x n )5(f)A(yi) ■ ■ ■ 54>A(y r , 







(2.7) 



(2i 



0A=O 



if the time component of anyone of the x's is larger than the time component of all the 
y's. Let us also remark that, by the definition of CTP generating functional, all time 
coordinates in the vertex functions are supposed to be positive so we can set, as well, 
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these functions to be zero for any negative time. Using eq. ()2.7|) together with eq. (J2.2j) 
one can write the equations of motions in the form 



5T 



(2.9) 



<M=o, <Pc=4> 



Notice that 2n-legs vertex functions with one a kg and In — 1 C legs are the only ones 
contributing to these equations of motion. Then eq. ()2.8|) guarantees that all the terms 
nonlocal in time in eq. ()2.9|) do satisfy causality. 

Perturbative calculations by diagrammatic expansion in the physical representation 
are based on the bare propagators 



G®(x,y) = G%\x,y) = -t6(y - x ) (*| 
G®(x,y) = G ( *\x,y) = -i9(x - y ) 



G { H A (x,y) = -2z(*|{0(y),0(x)}|*) c 



I free 



x), 4>{y)] l^) | free 

x)A(y)} 1^) \free 

, G<g(x,y) = 



(2.10) 



The bare retarded and advanced Green functions G^ and G^ do not depend on the 
initial state and are translational invariant. The bare correlation function G^ = jG AA: 
instead, does depend on |\&). The vertices are 




* ^0 Tijkn 



1 ^0 T ijkm 



i ^0 T ijkm v m 



= I ^0 Tiikm v n 



k j k j 

where solid lines represent C legs while dotted lines represent a legs. 

As stated above, the HF approximation consists in the resummation of diagrams with 
daisy and superdaisy topologies. By considering this diagrammatic resummation it is easy 
to verify that the corresponding effective action has the following general structure 

r H F[0A, 0c] = - (0a| □ |0 C ) - X, rj] 
where T is a functional of the following composite matrix fields 

iij{x) = Cii (x)0 cJ (x) , XiA x ) = 0c,i(a;)0A,i(a;) , Vij{x) 
We now introduce, for ease of notation, the new object 



0A,i(aO0A,jOr) 



1 bT 



2 5xtj(x) 
9 



x=v=o 



which is a functional of £ij(x) = 4> c ^(x)(f) C} j(x) = <pi(x)(pj(x) only. Then the equation of 
motion in the HF approximation takes the form 

= {D(Jy + 2^.[e]}0 i = O 

<t>A=0, 4>c=<t> 

We will consider this as a general form for mean-field-type background field equations 
and it will be that base for our definition of a modified HF approximation. 




3 The case N = 1 

We begin by considering a single scalar field theory with spontaneously broken Z2 sym- 
metry. The study of this simpler case provides valuable insight into the general features of 
the diagrammatic resummation performed by the HF approximation. In particular, this 
allows to understand the origin of the HF shortcomings with respect to renormalizability 
and RG-invariance and to determine the general recipe for the definition of a modified 
renormalizable and RG-invariant mean field approximation. Many of the results of this 
section will hold true also in the more general case of a theory with iV scalar fields. 



3.1 Analysis of the HF approximation 

For the N=l case the Hartree-Fock equations of motion [see eqs. (J2.6|) ] reduce to 
|n + m 2 + \\ Q i(x) + |A ^ ^Wj;)|j^) = 

\u p (x)\ 2 [> u k (x) = 



jn + mg + iAo£(a;) + !Ao 



I , -.12 



(2tt) 3 



'p 2 <A 

where we recall that = <ft 2 (x). Notice that, so far, these are still the equations of the 
regularized theory, written in terms of bare parameters (Ao, ttIq) and explicitly dependent 
on the sharp cut-off A. 

The phase with spontaneously broken symmetry, which is the subject of this paper, is 
defined by assuming the existence of a static homogeneous vacuum solution of eqs. (|3.1|) 
with nonzero background field <j>{x) = v. The corresponding mode functions, that will be 
named u^ vac \ have the following plane wave form 

u { r\x) = ^e^-^ , ul = e + m 2 (3.2) 

where m 2 = |A f 2 and the vacuum expectation value v of the background should satisfy 
the gap equation 

= ml + l\ v 2 + ^\ [ 7^^- (3.3) 
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The values of bare parameters for which this equation admits a nonzero solution v are 
those corresponding to spontaneously broken symmetry and are those we are considering 
here. 

We now introduce the mean field V according to 

V(a;)=mg + iAo« 2 + |AoAe(a;) + iAo / jA^ \u P (x)\ 2 (3.4) 

Jp 2 <A {^ 7T ) 

where A£ = £ — v 2 . Notice that V = on the static solution <p = v and u = u( va " c \ This 
definition allows us to rewrite eqs. (|3.1|) in the form 

{□- ±A o A£(x) + V(x)}0(x) = 

(3.5) 

{U + m 2 + V{x)}u k {x) = 

Recalling the definition of the background field equation in terms of the 1PI effective 
action 

^fk=o=( n + ^])<I> = Q (3-6) 

and comparing it with the first of eqs. ()3.5|) . we read out 

T%\ = \V-\\^i (3.7) 

Notice that V is regarded here as a functional of £ = <p 2 - The implicit dependence on £ is 
determined by solving the second equation in ()3.5j) with a generic background. To obtain 
a self-consistent equation for the functional V[£] from the HF equations of motion some 
preliminary definitions are required. 

• We introduce the free mode functions defined as solutions of the free equation 
with mass m 2 = | Ao v 2 and with the same initial conditions of the exact mode 
functions, that is 

(n + m 2 )u ( °\x) = , u { °\x,0) = u k (x,0) , ii ( £\x, 0) = u k (x, 0) 

• In terms of these we introduce the free correlation 



G^(x,y) = Re[ -^u$\x)v$\y) (3i 



Finally we define the free retarded and advanced Green functions as 



Gf{x-y) = Gf{y-x) = j 



d4 P 1 e -ip(x-y) 



(27r) 4 p 2 — m 2 + iepo 
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By using the above definitions we can cast the equation of motion for the mode functions 
into a convenient integral form 

u k {x)=u^\x) + J d 4 yG { ^{x-y)V{y)u k {y) 

In this equation, and everywhere else from now on, all field (0, mode functions, J 7 ', 
V, etc..) are to be thought as defined only for positive times (initial conditions are at 
the limit point t = + ) and all time integrations are restricted to positive values, as 
appropriate in an initial value problem. Moreover it is convenient to introduce a more 
compact operator notation with implicit space-time integration, that is 

u k = u k 0) + G ( ° ) Vu k (3.9) 

where the caret * turns a vector of the functional space into a multiplication operator 

" : v (x) — > v(x, y) = v(x) 5^\x — y) 

^From now on we shall use this notation throughout the paper. 
Eq. ()3.9|) can be formally solved for u k 

where 1 stands for the space time delta function 5^ (x — y). Then the cutoffed correlation 

d 3 p 

'|p|<a 
can be written as 

G = G[V] = [1 - G^V}- 1 G (0) [1 - VGf}- 1 (3.10) 

in terms of V, of the free retarded and advanced Green functions and of the free correlation 
function. In conclusion, by renaming the correlation at coincident points as 

I[V](x) =G[V](x,x) (3.11) 

and substituting into the definition of V in eq. ()3.4j) we obtain the sought self-consistent 
equation 

V = ml + l\ v 2 + ^\o(A^ + I[V]) (3.12) 

Notice that eq. ()3.12|) depends parametrically on the initial kernels (i.e. the mode func- 
tions initial conditions) through the explicit form of G^-°\ Now, before going any further in 
the discussion, we fix a particular choice for these kernels by considering the HF vacuum 
as initial state for the quantum fluctuations. That is to say that we start from equilibrium 
initial conditions for the mode functions 

u k (x, 0) = u k vac \x, 0) , u k (x, 0) = ii k vac \x, 0) 

which corresponds to the following choice of the initial kernels 

Q{k) = , S(k) = 

2u k 

By this choice some simplifying properties follow 



G{x, y) = Re / — — ^ u p (x) u p (y) 
J\p\<a i 2 V 
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• The free mode functions coincide with the vacuum mode functions w( vac ) at every 
time. 

• By the self-consistent equation we have V = at the point £(x) = v 2 

• The free correlation function defined in eq. (J3.8j) turns to be translationally invariant 
in spacetime. Explicitly 



Finally let us stress that, in spite of this choice, we are still considering an out-of- 
equilibrium problem since we allow for generic initial conditions for the background field. 
The use of different initial kernels will be discussed later. 

At this point some observations on the integral term I defined by eq. ()3.11|) and 
eq. (j3.1()|) are in order. 

1. / depends on the free retarded (and advanced) Green function and on the free corre- 
lation function. These propagators have the same form of their bare correspondents [see 
eq. (|2.1()jl ] but their mass is m 2 = ^XqV 2 rather than = tHq + ^XqV 2 . By using the gap 
eq. (J3.3)) . the relation between the two sets of propagators reads 



By recursively solving eq. 1)3.13)1 for the free propagators we can see that their definition 
corresponds to the resummation of all tadpole corrections. As we will see they play the 
role of internal dressed propagators in the diagrammatic resummation performed by the 
Hartree-Fock approximation. 

2. By expanding / in powers of V we can see that the term proportional to the n-th 
power contains loop with n + 1 free propagators. In particular the expansion up to the 
linear term reads 




= [1 - A A <4 0) /«] G?W [1 - | A Gf] 



(3.13) 




I[V](x)= G^ + G^VG^+G^VGf (x,x) + ... 



which, for later convenience, can be rewritten in the compact form 
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3. One can easily realize that diverges as A 2 and log A, diverges as log A and the 
loops with more than two propagators are convergent. It is therefore useful to introduce 
a new functional J containing only the convergent part of / according to 

j[V] = /[V] - J (1) - J (2) V (3.14) 



The previous observations lead us to rewrite the self-consistent eq. (J3.12)) in a particular 
"quasi-renormalized" form, more suitable to analyze the diagrammatic resummation of the 
effective action and determine its divergent graphs and subgraphs structure. The "quasi- 
renormalized" form is obtained by expanding the self-consistent eq. (|3.12j) around the 
point £(x) = v 2 (V = 0), explicitly solving the linear terms and writing a self-consistent 
equation for the higher A£ powers dependence of V. 

First we substitute eq. ()3.14|) into eq. ()3.12|) . applying the gap eq. ()3.3|) in order to 
simplify the constant terms. Then, by solving the terms proportional to V we obtain 



V=±£[A£ + J[V]] , 6 = \ [l-±A / (2) ] 1 (3.15) 
or, equivalently 

V=±#A£ + AV , AV= \6J[\6/\t + AV] (3.16) 

In conclusion we rewrite eq. (J3.7)) in terms of AV and of the explicit linear term as 

r = ±ft( 2 )A£ + ±AV , Q^ = f0-±A O 1 (3.17) 

We can now use eqs. (I3.1fij) and ()3.17j) to calculate the vertex functions relevant for the 
background equation of motion. In the physical representation of the CTP formalism 
these functions have the form 

with one 0a leg at the point x\ and n — 1 <p c legs at the points Xi (i> 1). 

First of all, from the general form of the background equation ()3.6|) we can see that 
the vertex functions (J3.18|) are built from the functional derivatives of T' w.r.t. £ at the 
point £ = v 2 . To ease notation we rename such variations by 



new..,*)-- S ™< 



5£(x 2 ) . ..5£(x k ) 

Some examples of the expressions of the vertex functions in terms of the are 
T {1 \x) =SlW(x)v , T ( - 2 \x ll x 2 ) = [□ + ft (1) (x 1 )]5 li2 + 2ftf 2 V 



123 



V ' 



( Xl , x 2 , x 3 ) = 2 fig 5 12 v + 2 fig 6 23 v + 2 fig 5 13 v + 4Q 

(3.19) 

r< 4 ) (xi, x 2 ,x 3 , x 4 ) = 2{ fig 5 12 5 U + fig 5 1A 5 23 + fig 5 13 5 2A } + 8 fig 34 v* + . . . 

■ ■ • + 4{ fig 4 5 12 + fig 4 5 13 + fig 3 s u + fig 4 s 23 + +fig 3 5 U + fig 3 5 2i } v 2 
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where we have introduced the shorthand notations 

= <$ (4) - Xk) , &fl k = ^ (2) (xi,..., x k ) 

/From eq. ()3.17|) we can see that = 0. So that = 0, which is just the statement 
that v is the vacuum solution of the background equation. Moreover, by eq. (j3.17|) we 
have 

= §0-§A Q l = §A Q l-fAo J< a > [1 - | A JW]- 1 (3.20) 

Expanding the above expression in powers of Ao we see that f^ 2 ^ is the sum of a classical 
term | Ao plus the resummation of all the chains of 1^ loops. 

The with k > 2 are obtained by using eq. ()3.17|) and the self-consistent equation 
for AV [eq. (j3.16|) ]. One can easily see that they are built up with loops having three or 
more propagators and are therefore finite in the limit A — > oo. The loops are attached 
to each other and to the external legs by the effective vertex 6 [see eq. (I3.15|) ]. Some 
graphical examples are 




The vertex 9, by its explicit definition in eq. ()3.15jl . is the sum of the classical term Ao 
plus chains of 1^ integrals 

^Aol-iA^l-iAo/^r 1 (3.21) 



3.2 Renormalization 

Let us now apply the standard renormalization procedure to the HF approximation of 
the effective action. As we will see it fails to render all the vertex functions finite. 

The first renormalization condition is simply the request that the v.e.v. v has a 
physical cutoff-independent value. Then the gap eq. ()3.3|) must be regarded as defining 
the bare mass parameter Wq as a function of Ao, v and the cut-off A 

mg(A ,t/,A) = -±A ot / 2 -±Ao/ J- (3.22) 

J P 2<A (27r) d 2u k 

Substituting the above function for the constant m\ in the equation of motion turns to 
be enough to remove all the A 2 dependence in the vertex functions. 

The second renormalization condition follows by fixing the value of the quartic coupling 
at some fixed energy scale. In order to do this we introduce the Fourier transform of the 
vertex functions [see eq. (|3.18j) ]. namely 



/m 
n=2 



(x 1} ... } x m )e i ^^ Pn - {xi - Xn) 
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where translational invariance, which is a consequence of our choice of initial conditions, 
has been assumed. For later convenience we define also the Fourier transform of 



J ^ — o 



We recall that, as appropriate to a causal initial value problem, all Fourier transforms are 
analytic in the upper complex po~halfplane. 

Now, in the present out-of-equilibrium context, we can define the symmetric point at 
the scale s by setting the momenta entering the C legs (p2, ps, p^ = —pi — P2 — P3) to 
the same, purely spatial, value 

~\pi =P2 =P3 =Pi = \q a , Qs = (0,q), \q\ = s (3-23) 

with an arbitrary direction q. The usual coupling renormalization condition is obtained 
by evaluating r^ 4 ^ at this point and requiring it to be equal to the renormalized coupling 
constant A 

A = f( 4 )(-§ 9s ,ig s ,ig s ) (3.24) 

Now, recalling the general expression of T^ 4 ^ in terms of the in eq. ()3.19|) , and 
performing the Fourier transform, we can rewrite eq. (|3.24|) as 

A = 6 OP>(q.) + 4{ ^\\q si \q s ) + 2 n®(-§fc, \q s ) + & 3 \-q s , ±g s ) 



+ n (3 \lq s , -q s )}v 2 + 8 n (4 \-lq s , |g s , \q s 



However, in our case we prefer to slightly change this standard procedure in order to avoid 
lengthy calculations and obtain a better comparison with the unbroken symmetry case. 
We therefore substitute eq. ()3.25j) simply with 

\ = QQ^(q s ) (3.26) 

We have omitted in this way the contributions to the quartic coupling originating from 
the superficially convergent and Hence the change from eq. (|3.25|) to eq. (J3.26)) 
is to be regarded as a legitimate finite part redefinition of the coupling constant. 

Now, by the HF expression for eq. ()3.20|) . the renormalized coupling constant 
explicitly reads 

1-|A„ /<«(«),) 

This has to be compared with the corresponding formula of unbroken symmetry case in 
ref. jnOl- Notice that in ref. [HUl we have adopted the opposite sign convention for I^ 2 \ 
Eq. (|3.27|) . once inverted, defines Ao as a function of A, v and (the logarithm of) A. 

As already stated above, the renormalization procedure just outlined fails to define a 
sensible renormalized theory. In particular we can individuate the following shortcomings. 

There is clearly a pathological behavior of the effective quartic coupling A as a function 
of the bare parameter A at fixed cut-off A. In fact A has the correct 1— loop Aq term 
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dictated by perturbation theory, but certainly fails at higher orders, since is exhibits an 
unphysical behaviour, growing to a maximum value at Ao = A™ ax , then decreasing to zero 
and to even more unphysical negative values (returning positive only for very large values 
of Ao). This implies the breakdown of the HF approximation for values of Ao greater than 
\max - n a th eor y a t fixed cut-off, as compared, for instance, with the standard 1— Loop- 
Renormalization-Group improved relation which reads 



1 o + ilogA + ... (3.28) 



A An 16tt 2 



where the dots stand for some suitable choice of the renormalization scale and of the finite 
parts. The relation in eq. (|3.28|) is monotonically increasing with Ao, at fixed A, to an 
asymptotic plateau value X max . As a consequence it can be inverted determining a single 
bare coupling Ao value for any A and A < X max . A full trajectory of Ao as a function of 
A, up to the vertical asymptote at the Landau Pole value A LP , corresponds to a single 
renormalized theory describing the same dynamics at momentum scales much smaller than 
A. On the other hand, in the Hartree-Fock approximation, the non-monotonic behaviour 
of the relation in eq. (|3.27|) spoils this one-to-one correspondence between bare and 
renormalized parameters (at fixed cut-off) which holds true only for small coupling. 

Moreover, we see that imposing a finite value to A at the chosen scale s fails to render 
finite the running coupling (i.e. the generalization of eq. ()3.26|) to any value of momentum) 

A(p) = 6Q (2 \p) (3.29) 
at any momentum p ^ q s and even at any q s i with s' ^ s. 

For what concerns the higher order terms f^ fc ) with k > 2 we can see that a logarithmic 
dependence on the cut-off persists even after imposing the renormalization condition in 
eq. (|3.26|) . That is they do not parametrically depend solely on v and A, but also on A 
and therefore on log A. This happens for two specific reasons 

1. The internal propagators G^ , G^ and G^ introduce an explicit and non removable 
dependence on Ao due to the presence of the resummed mass m 2 = | Aof 2 . Notice 
that, unlike the unbroken symmetry case, here the resummed contribution to the 
internal propagators differs from those to the external propagator defined as the 
functional inverse of the two-legs vertex function. To the latter contribute also 
diagrams as 




as can be seen from its explicit form [see eq. (j3.19|) ] 

f ( 2 \p) = -p 2 + \ A(p) v 2 = -pi + m *-l \l p)(p) [l - i Ao I {2 \p)]- 1 v 2 (3.30) 

This is finite at momentum q s while, according with what we said above about X(p), 
logarithmic divergences appear as p ^ q s . 
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2. In the HF definition of with k > 2 there appears the effective vertex 9. By the 
relation [see eq. ()3.2H) ] 

6{p) = | A(p) + | A 

we can see that, after imposing the renormalization condition, an unresolved Ao de- 
pendence persists even when p = q s . For p ^ q s further cut-off dependences appears 
due the problems concerning the renormalized running coupling in eq. (|3.29j) . 

We conclude that the the HF approximation cannot be renormalized by the standard 
renormalization procedure: there is an unphysical bare-to-renormalized coupling relation 
and a plain failure to eliminate divergences in the subgraphs of the resummation. In 
the next subsection we will define a modified HF resummation by explicitly requiring 
renormalizability and a 1— loop-renormalization-group improved relation between A e 
Ao [see eq. (|3.28|) ]. By comparing the two resummations we will recognize the cause 
of the renormalization problems of the Hartree-Fock approximation in the incomplete 
resummation of Leading Logarithms of the cut-off. 

3.3 Improved HF approximation 

Our recipe for the improvement of the HF resummation consists in two fundamental steps 

1. Fix some features, proper of the HF approximation, that we want to maintain 
throughout the modification since they provide a minimal definition of a mean field 
resummation. The results of sec. 13.11 provide a parameterization of the class of 
approximations having such defining features. 

2. Require explicitly renormalizability and RG-invariance in order to fix the form of 
the arbitrary parameters. 

In conclusion, once the result is obtained, we will be able to establish a diagrammatic 
interpretation of our modified approximation. 

The main features of the HF approximation are encoded in 

• The general mean field form of the background equation of motion 

5Yi m 



6a=0 



(n + 2F'[Z](x))(t)(x) = (3.31; 



The self consistent definition of T' (through the mean field V) 

T = [|fi (2) -±#]A£+§V , V = \6A£ + \6J[V} (3.32) 
where we recall the definition of J by eq. ()3.14|) and eq. ()3.11|) . 
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As we have seen in sec. 13.11 the relations in eqs. ()3.31|) and ()3.32|) provide a general 
recipe for building up all vertex functions using f^ 2 ) , the effective vertex 9 and the free 
propagators. Moreover, the self-consistent eqs. ()3.32|) for T' correspond to the following 
general form of the equation of motion for the mode functions 

{n + m 2 + V(x)}u k (x) = , V = i#{A£ + /-/ {1) -I&V} 

with 

i(x)= L„w? Kix)? 

We now regard QS 2 \ 9 and the effective mass m 2 of the free propagators as tunable 
parameters defining a class of approximations that share the same general diagrammatic 
structure. In other words, we abandon the HF definitions for these parameters, that read 

m 2 = \\ v 2 , 9 = Ac [1-iAo/^f 1 , ft (2) = ±fl-§A l (3.33) 

and look instead for new definitions that ensure proper renormalizability and RG-invariance. 

We can actually further specify the form of the tunable parameters. By looking at the 
HF definitions in eq. (J3.33)) of Vt^ and 9 we see that they must have the general leading 
log structure 

^^AoiMAo/^) , fl = A F 2 (A / (2) ) (3.34) 

in terms of two functions, F\ and F 2 , of a single variable (the evaluation of these functions 
on the operator is obvious if we consider the Fourier transform). We assume that 
the same holds true for m 2 . This can be understood looking at the explicit form of the 
two-legs function in eq. (|3.3(J|) . We consider m 2 as a function of evaluated at zero 
momentum since we don't want to introduce any new mass scale dependence, that is 

m 2 = A F 3 (Ao/ {2) (0))*; 2 (3.35) 

Notice that changing m 2 corresponds to changing the equilibrium solution according to 
eq. 0S3). 

To conclude we observe that the HF approximation resums correctly the 1— loop per- 
turbative order. To maintain this feature we should require that Vt^ matches at tree 
and 1— loop order while 9 and m 2 match at tree level. Explicitly we have the following 
conditions on the F's 

= \ + \x + 0(x 2 ) , F 2 (x) = 1 + 0(x) , F 3 (x) = \ + 0{x) , (3.36) 



Now we have to define our modified HF approximation by fixing new explicit definitions 
for the functions F's. We obtain renormalizability with the correct 1— loop beta function, 
by assuming the following logarithmic dependence for the bare coupling Ao 

dX ° 3 r A 2 + 0(A- 1 ) 



d log A 16?r 2 
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and requiring that the parameters ^ 2 \ 9 and m 2 do not depend on log A. Hence, differ- 
entiating eqs. (J3.34J) and ()3.35|) with respect to log A, yields 

(l-|x)F J f(a;) + |F J (a;) = 0->F J (x) = ^-^- 7=1,2,3 (3.37) 
where Aj are integration constants and we have used 

ai3iA = -8^ 1 + °( A " l » 

Eqs. f!3.37|) have a unique solution that fulfill the matching constraints in eqs. ()3.36|) . 

F 2 {x) = 6F 1 (x) = 3F 3 (x) = —L- 

In terms of the parameters Q( 2 \ 9, m 2 

e = Q^ = X [1-I Xol^}- 1 , m 2 = l\ [l-IX I {2 \0)}- 1 v 2 (3.38) 
These are the explicit forms of the parameters in our improved HF approximation. 

Now applying the renormalization condition in eq. ()3.2fij) with the new definition of f^ 2 ) 
we obtain the following bare-to-renormalized relation 

A = 5-^= (3.39) 

l-IX F 2 )(q s ) 

Comparing this with eq. ()3.28|) we can see that our procedure has reproduced the correct 
1— loop-renormalization-group improved behaviour while at the same time fixing the 
finite parts. Substituting eq. (|3.39|) into eqs. (|3.34|) and eq. (|3.34j) we obtain a manifestly 
finite form for the parameters 

e = 6£lW=\[l-l\J®]- 1 , m 2 = |A[l-|AJ( 2 )(0r 1 t; 2 (3.40) 

in terms of A and of the subtracted integral 

j(2)=/(2)_/(2)(g a )l 

In particular, recalling the natural definition of the running coupling constant in eq. (|3.29|) . 
we see that now 

6(p) = A(p) = - 3 ^~. 2 , . . (3.41) 
1 - I A J( 2 ) (p) 

Now, for what concerns RG-invariance, one verifies from eq. ()3.39j) that the parameteri- 
zation of Ao does not depend on the renormalization scale. In fact, renormalizing at scale 
s with constant A = X(q s ) or at scale s' with constant A' = X(q s r) indeed defines the same 
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Therefore eqs. (J3.38|) can be regarded as manifestly RG-invariant definitions of the pa- 
rameters. 

Finally let us observe that, since depends on m 2 , we have here an implicit definition 
of the physical mass by the finite self-consistent relation in eq. (|3.4U|) . rather than an 
explicit definition, such as the tree-level m 2 = |At> 2 , which can be recovered only if 
s = 0. 

We give now a brief diagrammatic interpretation of the results just derived. First of 
all, with the new form of m 2 in eq. 1)3.38)1 . the free propagators are no longer defined as 
the resummation of tadpole corrections. In fact are included contributions from graphs 
as in fig.. These were actually present in the HF definition of the external propagator. 
Moreover contributions from graphs as 




are included in such a way that the log A dependence of m 2 (at fixed Ao) corresponds 
to the correct 1— loop-renormalization-group improved series. Notice that corrections of 
this type, if included completely, would define a momentum dependent self-energy. On 
the other hand it is implicit in a mean-field approximation that the internal propagators 
have a momentum-independent self-energy. Thus the contributions of these diagrams are 
to be included in a "local" fashion (i.e. with no momentum passing through the loops). 

Now, for what concerns the new form of in eq. (|3.38j) . we can see that it corre- 
sponds to the inclusion, in addition to the chain diagrams of the pure HF approximation, 
of all the Leading Log contributions from diagrams as 




+ sym + . . . 



As before their log A dependence is taken while their finite part is fixed by our procedure 
to be the same of 

In conclusion, the modification of the effective vertex 6 in eq. (j3.38j) corresponds 
to include, in the diagrammatic resummation that defines the fiw (k > 2), Leading 
Logarithmic contributions from graphs of the form 

The finite parts of these diagrams are chosen in such a way to maintain the main Hartree- 
Fock-like features of the effective vertex, namely the single channel structure and its form 
as a function of I^ 2 \ 
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Let us make some final observations on the obtained result. The equations of motion of 
this modified Hartree-Fock approximation read 

{D + V(x) -±#A£(x)}0(x) = , {D + \9(0)v 2 + V(x)}u k (x) = 

' (3.42) 

V = \ 6 { A£ + 1 — - J (2) V} 

As we can see, in order to obtain renormalizability and RG-invariance, we had to introduce 
space time non-locality in the equations of motion. Causality in this nonlocal evolution 
is guaranteed by the analyticity in the upper halfplane of I^(p). Actually in the third 
equation in eqs. (|3.42|) the definition of V is implicit and it should be solved for V in order 
to have a manifestly causal form. This can be done paying the price of losing manifest 
finiteness 

V = \ A [1 - A (J (2) - f J (2) (g s ) l)]- 1 {A£ + / - /W} 



3.4 Other initial states 

We already pointed out, in sec. 13.11 that eqs. ([3.42 1) . with vacuum initial kernels and 
generic initial conditions on the background field, already describe an out-of-equilibrium 
problem. In this section we study whether and how we can choose different initial condi- 
tions for the mode functions without spoiling the properties of renormalizability and RG 
invariance. We already considered this problem in [3U| treating the unbroken symmetry 
case. The results in this section will be very similar. 

As an example consider an initial state of the same form of the HF vacuum 



G{k) 



2Qi 



S(k) = 



k 2 + M 2 



but with a mass M different from the equilibrium mass m defined in eq. (J3.38|) with 
M m. The corresponding free mode functions have no longer the plane wave form 



u (0) = ^(vac) [ gee eq _ 



U 



(0), 



,,ikx 



x) 



2 V / 2T^ 
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1 + — 


e ~iuJ k t _|_ 
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e ~iu> k t !_ 




0J k . 









and the free correlation function is not translationally invariant in time. Then the integral 
term in eqs. ()3.42|) 

7_/(D_j(2) V 

in spite of the subtractions still contains the superficially divergent contribution 



d 3 k 
(27P 



m 



M 2 



COS 2uJy.t 



that indeed diverges with the cut-off when t = 0. These initial time singularities can 
be removed by a Bogoliubov transformation on the initial state as first observed in 
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This transformation redefines the initial kernel in such a way that the leading terms of 
an high-momentum expansion are the same as at equilibrium 



1 tti 2 
Q(k) ~ — = + , + . . . 

2V¥ 4(fc 2 ) 3 / 2 

Initial singularities as well as any other divergence turn to be absent for any choice of kernel 
having the above large k expansion. A simple interpretation is that the renormalization 
procedure ensures finiteness for any initial Gaussian state belonging to the same Fock 
space of the HF vacuum. Extrapolating from this simple example a generic condition on 
the short-distance behaviour of the initial state kernel, we have 

1 m 2 
G(x, V) ^ -r-2j r + ^ lo S \ x ~ V + • • • 

This ensures the cancellation of all divergent terms are guaranteed by mass and coupling 
constant renormalization. 



4 The case N > 1 

We are now ready to consider the more general case of a scalar field theory with sponta- 
neously broken O(N) symmetry. In doing this we proceed following closely sec. 

Before we begin we fix some notational conventions that will help us to handle O(N) 
index structures while keeping formulas simple and similar to those in sec. El In particular 
we will use the standard matrix notation for objects with one and two indices 



[MV] l = M l3 V 3 , [M« M< 2 >] y = M® M$ 



(4.1) 

When objects with four indices are concerned, we will use the following conventions 
[T M] tJ = T ljkm M km , [T« T% jkm = T$ s T® m 

(4-2) 

tlijkm — 2\3ikfijm + &im$jk) > T T = TT =11 

where the last relation is regarded as restricted to tensors with the symmetry Tij km = 
Tijmk = Tjikm- Moreover we will use also the following definitions 

l(x,y) = 15^(x-y) , Il(x, y) = U 5 {4 \x - y) 
4.1 Analysis of the HF approximation 

The Hartree-Fock equations of motion for a scalar O(N) theory are shown in eqs. (J2.fij) . 
As a first example of the use of the notational conventions in eq. (j4.1|) and eq. (|4.2j) we 
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rewrite them in a compact indiceless form 



j[n + mo + §Aotr£]l + ±A T u p u p (x)^ </> 

|[D + mo]l + |A Q T £ + J 2 u p u p | u k = 



(4.3) 



where we recall that = 4>(x) <j) T (x) and the definition of r in eq. ()2.5|) . 

Exactly as in sec. 13.11 we select the values of bare parameters corresponding to the 
broken symmetry phase by requiring the existence of a vacuum solution with non zero 
constant and uniform background field, <f>(x) = v. Notice that here v is a O(N). The 
direction v = v/Vv 2 of the vacuum background field provides the definition of longitudinal 
and transverse projectors 

P L = vv T , P T = 1-P L 

The vacuum mode functions now are 

u k \ x ) = (2^fc) -1 ^ 2 exp[z(fc • x) 1 — iuJkt] 

= [{2u ktL )- 1/2 e~ iuJk ' Lt P L + {2u ktT y 1 / 2 e - luJk ' Tt P T ] e ik x 

where in the first line the power and the exponentiation are operations on matrices and 
u), u>l and ujt are defined as 

^l = ^l,L P L + ^l,T P T = k 2 l + M 2 , M 2 = m 2 L P L + m 2 T P T 

The longitudinal mass has the same value as in the N = 1 case, namely m 2 L = ^Xov 2 , 
while the transverse squared mass m^, as a function of Ao, v 2 and A, is obtained by solving 
the self-consistent equation 

= m 2 T + l\ [ -^-A— — ) (4.4) 

The vacuum expectation value v satisfies a gap equation which generalizes the one in 
eq. O 

V = ml + \\ Q v 2 + \\J *PJ— + ftN-l)* [ ( 45 ) 

V<A OO 2uJ k,L b Jp2<A (2tt) 3 2w fciT 

Notice that the gap equation involves only v 2 . In fact varying the direction of v we obtain 
different static solutions that correspond to the equivalent distinct vacua of the theory in 
the broken symmetry phase. 

The mean field V [see eq. (|3.4j) ] is now an object with two O(N) indices, defined by 
V(x) = m 2 1 - M 2 + l\ r[ax) + J ^ ^ \u p (x)\ 2 ] 
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Notice that again V = on the vacuum solution. In terms of V eqs. (|4.3|) read 

{□1 + M 2 -§A O £ + V}0(x) = 

{□1 + M 2 + V}u k (x) =0 

Recalling the general mean field form of the background field equation [T' here is a O(N) 
matrix) 

= (□1 + 2^])0 = O (4.6) 

we have 

r[^] = |V + |M 2 -|A e 
which should be compared with eq. ()3.7|) in sec. 13.11 

In order to derive a self-consistent equation for V like the one in eq. (|3.12|) . as in 
sec. 13.11 we introduce the free mode functions 

(□ + M 2 )uf (x) = , 4 0) («, °) = 0) > 4 0) («, °) = Wfc(a, 0) 
and define the free propagators 

G^ - y) = - y)P L + G% T (x - y)P T = / ^ ^gg^ 

G%(x-y) = G ( g ji (y-x) (4-7) 
(*,!/) = G? ) (x I y)P L , y + 4 0) (x I y)i? r ,« = Re / ^u%{x)u%{y) 

J\p\<A ( 27r ) 

Performing the same manipulations on the mode function equations as in sec. 13. II we can 
rewrite the cutoffed correlation, that now reads 

Gij(x, y) = Re / — — u pM (x) u P , a j(y) 
J\p\<a ( 27r ) 

as a functional of V 

G = G[V\ = [i - G^V]- 1 G (0) [i - VG^]" 1 
In conclusion we have the self-consistent equation for V 

V = m 2 l-M 2 + iA r(£ + /[V]) , I[V](x) = G[V](x, x) (4.8) 
This should be compared with eq. (|3.11|) and eq. (|3.12jl 

We keep on following closely the sec. 13.11 by fixing vacuum initial conditions for the 
mode functions. By this choice follows that for all times and the free corre- 

lation is translationally invariant. 
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Moreover we can repeat, with little changes, the observations made in sec. 13.11 about 
the structure of the integral I[V]. We can still interpret the free Green functions as 
effective internal propagators obtained by resumming all tadpole corrections to the bare 
ones. We can introduce the tadpole integral 

J« = G«>\x,x)=l£ ) P L + l£ ) P T 

and the two propagators loop integral 

J$L(* -v) = g rU* - y) y) + cgU* - y ) G<g(x, y ) = 

= if I {x - y) P L , ik P LJm + I®{x - y) P T , ik P T , 3 m+ (4-9) 



+ \ ItI( X ~ y) ( PL,ik PT,jm + Pr,ik Pl 



In terms of these we can define a functional J by subtracting to I its constant part and the 
term linear in V. By the same arguments used in sec. 13. II we can see that J contains only 
the cut-off convergent part of I. The definition of the J functional is formally unchanged 
respect to eq. (j3.14|) [here O(N) indices contraction according to the rules in eq. (|4.2jl is 
understood] 

J[V] = /[V] - J (1) - J (2) V (4.10) 

Notice also that, by using the two propagator loop integral of eq. ()4.9jl . we can cast the 
self-consistent definition of m\ in eq. f)4.4j) into a more compact form 

m\ = -|Ao4 2 2(0)m|[l - i Aq/SCO)]" 1 (4.11) 
The "quasi-renormalized" form is obtained by manipulations identical to those in sec. 13.11 

V = ±#A£ + AV , AV= \9 J[AV + ±0A£] , T' = \ + \ Q {2) A^ + \ AV 
The matrix f^ 1 ) and the four-indices objects and 9 are defined as 

nW=m 2 T P T , 9 = A [ft- iA / (2) j _1 , ft (2) = |Aofl (4.12) 

The vertex functions contributing to the equation of motion [see eq. ()3.18j) ] are ob- 
tained from the variation of eq. (|4.fij) with respect to C at = v. They can be expressed 
in terms of functional derivatives of T' 



as in eq. f)3.19j) . Notice that VL^ is completely transverse [see eq. (|4.12|) ] which implies 
that T^ 1 ^ = Vt^v = 0, namely the statement that the vector v is the vacuum static 
solution. 

Of course all the diagrammatic interpretations made in sec. 13.11 still hold true, now 
with diagrams carrying O(N) indices. 
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Before going further, we should introduce some notations that will be useful later 
on. In fact we will have to deal with four-indices tensors which have some particular 
symmetry properties and are functions of v (and of other scalar quantities) 

T^jkmsP) Tjj i l~ m (v) ^~ij'mfc(^') -^fcrmjl^) 

as, for example, 9 and tt( 2 \ Such tensorial objects admit a general decomposition 

T ijkm (v) = T Q (v 2 ) t? jkm (v) a = 1, . . . , 5 (4.13) 
in terms of the five elementary tensors 

tijkm = ViVjVkV m , tijkm = \ {Pr,ik PT,jm + ^T,im ^Tj'fc) 

tijkm = PTjj Pr,km , t ijkm = {Ph,ij Pr,km + Pt,%j PL,km) 

tijkm = \ {Pl,ik PT,jm + Ph,im Pr,jk + Pr,ik Ph,jm + PT,im Ph,jk) 

Notice that the coefficients of the decomposition are functions of v 2 alone. As an example, 
by eq. (|4.12)) the coefficients 9 a of the decomposition of 9 are 

6 X = A [1 - l(N + 2) A J#] {1 - i(JV + 1) A - |A [1 - |(JV + 2) A ig] I" 1 

02 = |Ao [1 — | Ao 4j.] 1 , #4 = |Ao [1 — I Ao 1 

5 = 1 - |(iV + 2) Ao iS]- 1 , 3 = 05 [1 - § Ao l| 2 2] [i - | Ao ig]- 1 

(4.14) 

Notice that using eqs. (|4.14jl the compact form of the self-consistent definition of m 2 ^ in 
eq. 1)4.11)1 can be rewritten as 

m% = ±\ v 2 - 1 9^(0) v 2 (4.15) 

The coefficient of Q can be determined from those of 9 by eq. ()4.12|) and the decomposition 
of the II tensor 

4.2 Renormalization 

Now we apply the standard renormalization procedure as we did in sec. 13. 21 for the N = 1 
case. We will see that the same problems are present also in this case and some others 
will arise due the presence of the spontaneously broken continuous O(N) symmetry. 

By fixing the equilibrium value of the background field we provide the first renormal- 
ization condition that defines the bare mass as a function of Ao, v and A 

m 2 (Ao,,;,A) = 4W-!Ao / ^ -J_ + I(JV-l)Ao / i% 

V<a OO 2uJ k,L b y p 2 <A (2tt) 3 2a; fc ,T 



g 9 a — g Aq l 



a 



1,2,4 



\9 r 



a 



3,5 



(4.16) 
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We recall that m| = |Aof 2 and is given (in function of Ao, v and A) as solution 
eq. (|4.4jl . Again, enforcing this condition is enough to remove all A 2 dependence from 
the vertex functions. 

The second renormalization condition is conventionally obtained by evaluating at the 
symmetric point [see eq. ([3.23)1 ] the "all-longitudinal" component of the four-legs vertex 
function r^ 4 ^ and requiring it to be equal to the renormalized coupling A 

A = f§L(-f Qs, \ q„ \ q,) h Vj v k v m 

Explicitly expressed in terms of the variations of T' the above expression involves the 
(longitudinal components) of f^ 2 **, f^ 3 -* and f^ 4 -* . As in sec. 13.21 [see eq. (I3.2fij) ] we drop 
the contributions from Vl^ and f^ 4 ** obtaining a simpler definition 

* = 6 ^kjQs) Vi Vj v k v m = 6 nf\q s ) (4.17) 

or, thanks to eq. (j4.12|) . 

A = 3 A [1 - ^ A Itt(Qs)} i < A lg v 

This has to be compared with the corresponding relation of the N — 1 case in eq. ()3.27j1 . 

First of all eq. ()4.18|) shows the same pathological dependence of A on A , at fixed A 
and v, that we found in sec. 13.21 As already discussed this prevents a consistent map 
between bare parameters and renormalized ones unless we restrict to a small coupling 
regime. The correct 1-Loop-Renormalization-Group improved relation, instead, reads 

where again the dots stand for some suitable choice of the renormalization scale and of 
the finite parts. 

For what concerns the cut-off dependence of f^ 2 ** we can see that things are even 
worse than in the N — 1 case. In fact, not only the log A dependence is not removed by 
the renormalization condition from (l^(p) for p ^ q s , but also, for components fl a with 
a ^ 1 logarithmic divergences already appear at p = q s as can be verified by the explicit 
definitions in eqs. ()4.16j) . 

Regarding the cut-off dependence of the VL^ with k > 2 the same observations of 
sec. 13.21 still hold in this case. In particular we still have log A dependence in the masses 
of the free propagators. This is immediate for m 2 L = | Ao v 2 and can be easily verified for 
m\ by considering its self-consistent definition in eq. (|4.4j) . Moreover, as in the N = 1 
case, log A divergences persist in the effective vertex 9 after imposing the renormalization 
condition [eq. (j4.17j) ]. 

Besides these renormalization problem analogous to the N — 1 case, another very 
important aspect, which is peculiar of the > 1 case, must be pointed out. First of all, 
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we may verify the formal O(N) invariance of the HF approximation also in the current 
out-of-equilibrium context. In fact, eqs. ()4.3|) are manifestly 0(N) symmetric and so it 
should be for the effective action 

r 1P i[fa,<l> c ]=r 1 p I [R<l> A ,R<l> c ] , VReO(N) 

Hence for infinitesimal R— 1 + e— 1 — e T , to first order we must have 



= / d 4 x <{ 4>l{x) e ( — r 1P7 [0 A , <f> c ) ) + f c (x) e ( -^j- V 



1PI\<PA, 



!>a(z) / \5(/) c (x) 

Upon variation with respect to <p A at (fi A = 0, <p c = v and recalling that I^ 1 ) = 0, we 
conclude 

J d A xv i e ij T i $(x-y)=0 f^ 2) (0) = 

This Ward identity, stating the masslessness of the external transverse propagator (the 
inverse of T^\p)), is indeed satisfied by the HF approximation 

r( 2 ) = rf> p L + vf p T 

tf(p) = - P 2 + 2^\ P )v 2 (4.20) 



fP(p) = -p 2 + 0« + nf\p) v 2 = -p 2 + m 2 T + \ h(p) - \ Ac 



v 2 



thanks to the self-consistent definition of the transverse mass, eq. ([4. 15)1 . But the trans- 
verse mass itself, which enters the internal transverse propagators G$ T and Gj) does 
not vanish at all (it actually diverges as log A), preventing a consistent interpretation 
of the transverse modes as Goldstone bosons. This is a well known problem of the HF 
approximation (see for instance ref. |32j). 



4.3 A class of improved HF approximations 

We shall now try and improve the HF approximation to recover the correct properties of 
Renormalizability and RG-invariance. We shall also require that this improved resum- 
mation is gapless; that is, we shall impose that the internal transverse propagators are 
massless. To this end we follow closely the procedure of sec. 13.31 while stressing some 
important new features which appear due to the continuous O(N) symmetry and the 
presence of two kinds of fields (transverse and longitudinal) with different masses. Let us 
anticipate that in this case our procedure does not single out a unique solution but rather 
an extended class of resummations that share all the required properties. 

As in sec. 13.31 the first step consists in fixing some fundamental properties of the HF 
approximation that we want to preserve. First of all the general structure of the diagrams 
resummation which is encoded in the general mean field form of the background equations 
of motion: 

= (Dl + 2f'i^W = (4.21) 
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and in the self-consistent definition of T' 

r=!fi (1) + [ift (2) -i#]A£+iV , V =\6^ + \6J{V] (4.22) 

which imply the following mean-field-type equations of motion for the mode functions 

[□1 + M 2 + V]u k = , V= ±0{A£ + /-/ (1) -J (2) V} (4.23) 

As in sec. 13. 31 we now regard fiW, fiW, 6 1 and the free propagators masses m 2 and m 2 - (in 
matrix form M 2 ) as tunable parameters that we are going to change w.r.t. to their HF 
definitions 

e = \ \n-i\ m- 1 , n^=m 2 T p T , = \e- ±a tt 

3 (4.24) 

m 2 = lA oW 2 , m 2 = 1\ qV 2 - \6 A {®)v 2 

The first important difference with the N — 1 case is that now f^ 2 ', and M are not 
independent, but must fulfill certain O(N) symmetry constraints. In fact the form of 
the equations of motion in eqs. (|4.21|) and (|4.23|) is not manifestly O(N) symmetric for 
generic tunable parameters. More precisely they are covariant under contemporaneous 
rotations of the field <p, of the mode functions u p and of the vacuum expectation value 
v. The latter, in fact, enters in the definitions the propagators (see the decomposition in 
eq. (|4.7)l ) and of VL^ and (see the general decomposition rules in eq. (|4.13jl ). We have to 
explicitly require symmetry under rotations of <fi and u p alone or, equivalently, invariance 
under rotations of v. In other words we must impose the proper O(N) Ward identities. 

Let us therefore define the operator generating the infinitesimal variation under an 
O(N) rotation of v 



where e is an arbitrary antisymmetric matrix. It is immediate to verify that it annihilate 
any invariant function. 

First of all, to require invariance in the mode functions equations (|4.23|) . we have to 
impose 

5[M 2 + V] = (4.25) 

Solving explicitly for V the second equation in eqs. ()4.23|) we can rewrite the above 
condition as 

= SM 2 + S{ [1 + 10 1®}- 1 [Af + 1 - (4.26) 

Notice that £ and Uk are independent and by hypothesis invariant variables (i.e. they 
must depend only on v 2 ). Then eq. (|4.2fi|) implies 

5{[1 + \6 I^Y 1 6} = Q , = 5M 2 + [l + \ei m ]- l e[-v 2 8P L - 51 {l) }} 

By some easy algebraic manipulations we can rewrite the first equation above as 

56= \ 6 5I {2) 6 
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which has a simple and natural diagrammatic interpretation. Making the tensorial struc- 
ture explicit, we obtain the following conditions on the coefficients 6 a 

6» 2 = 6> 4 [1 - \ (i^rp - ) ] _1 

e 5 = [01-04+1 *i ( J S - 42) ] [i + 1 ^4 ( 4? - 41 ) r 1 ( 4 - 27 ) 
03 = 05 [i + 104(42 - 4 2 l)][i + 10 4 (4? - 4 2 !)]- 1 

Notice that these relations uniquely fix the form of once 9\ and 84 have been specified. 
Now, using these relations and the following rules for the explicit variations 

5M 2 = (m 2 L - m 2 T ) 5P L , 51® = - J^) 5P L = (m| - m|) 42(0) <5P L 

we can rewrite the second equation in eqs. ()4.3|) as 

m 2 L -m 2 T = ±v 2 4 (O) (4.28) 

In the second step, we impose invariance of the background field equation (j4.21|) by 
requiring 

5T' = (4.29) 
This can be written, by using eqs. (J4.22|) and (J4.25j) and the independence of £ and u, as 

5^ = \56 , 5Q« -5M 2 -(^ -\6)v 2 5P L 

Explicit tensorial calculation for the first equation yields 

f2 2 2 = £l± + \ [02 — 04] 

= - + 1[0 3 - e x + e 4 ] (4.30) 
nf ] = nf ] + l[0 5 -0i + 4 ] 

While for the second we obtain 

= -nf\o)v 2 p T (4.31) 

In conclusion, according to eqs. ()4.27|) . ()4.28j) . ()4.3()|1 and ()4.31|1 . the actual independent 
parameters are 0i, 4 , m|, and fi^ . The integrated versions of eqs. (|4.25j) and (|4.29|) 
are obtained by writing M 2 + V and T' in a manifestly O(N) symmetric form 

m 2 + v = pWi + | 7 (1) [£ + A , r = ^ 1 + i 7 {2) e + 1 [M 2 + V] 

7i$L = 71^ Sijhm + \ li^ [dikdjm + Sim$jk] , a =1,2 
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with 7 < - Q ) and denned in terms of the free parameters as follows 

= [01-9,-1 e, {ifl - 4 2 l)] [i + § 4 (4 2 2 - 4 2 })} [i _ i <? 4 4 2 )]-i a- 1 



2 

,(i) _ a n _l i r( 2 )i-i 



(4.32) 



a = [i + 1 4 [ /S - ni$]] [1 + 19, if L ] -a=±ei 4?[i + §04 4 2 2r 1 

7 « = 4 [i + i4 2 l] 

^ = m| - |#(0) b 2 + /«] - !#(()) [u 2 + /« + (AT - 1) jW] 

One can verify that the standard HF definitions of the tunable parameters indeed sat- 
isfy eqs. (I4.27|) . (|4.28|) . ()4.3U|) and (|4.31j) and that the corresponding parameters of the 
manifestly O(N) symmetric form are 

7 « = 2 7 «=4 7 f = |A 1, 7 ! 2) = 0, ^ = \ml, ^ = 

Now, as in sec. 13. 3| we require that 9 and f^ 2 ** have the same general leading log structure 
characteristic of their HF definitions [see eqs. (|4.14jl and eqs. (j4.1fij) ]. namely 

^i 2 = Ao-Fi(Ao4l' ^o4l' Ao4t! N) , 9i = A F 2 (Ao4l' <Wtl> Ao4t! N) ^ 
= Ao-^3(Ao4l> Ao4l> Ao4tj N) , 6> 4 = X F 4 (Ao4l! Ao4l> Ao4t> -A") 

where the y, z; AT) are generic functions of commuting arguments since the operators 

4^) 4l an d Itt are diagonal in Fourier space. 

Notice that the same structure is inherited by all the others components of 9 and fl^ 
as can be verified by eqs. ()4.27|) and eqs. ([4.30)1 . We assume a similar structure also for 
m\ as we did in eq. ()3.35|) for the N — 1 case 

m 2 L = A ^ 2 F 5 (Ao4 2 i(0), A jg(0), A o ig(0); iV) (4.34) 

Notice that, by means of eq. (|4.28|) . the same structure is inherited by m\ which is 
consistent with its HF definition [see eqs. ()4.24|) ]. 

We have thus reduced our parameterization freedom to five functions Fa of three 
variables. We still have to require that these functions fulfill some further properties that 
hold true in the HF approximation. 

Hartree-Fock matches perturbatively at 1-loop order. This matching requires that 9 
and ml match at tree level and that matches at tree level and 1-loop order. Actually 

(2) 

Q\ must match only at 1-loop level since the tree level term never appears in the vertex 
functions. Explicitly we have the following conditions on the F's 

F 1 (x,y,z;N) = ± + \x+±(N-l)z + ... , F 2 (x, y, z; N) = 1 + . . . 
F 3 (x,y,z;N) = c + \y + ... , F 4 (z, y, z; N) = § + . . . 

F 5 {x,y,z;N) = § + ... 
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where Co is a purely numerical arbitrary constant. One can easily check that the above 
conditions, together with the symmetry relations in eqs. (|4.27|) . ([4.28)1 . ()4.30|) and (|4.31|) . 
are enough to guarantee the matching at 1-loop order of all the vertex functions. 

In the N — > oo limit the HF resummation reproduces correctly the equations of the 
usual large N approximation. This can be done as follows. First we restrict to the special 
case of a background field which maintains a fixed direction (i.e. the direction of the 
vacuum expectation value v) , that is (fii(x) = (p(x) -Oj. Then we can reduce the equations 
of motion into a projected form by setting 

u k = u kL P L + u kT P T , V = V L P L + V T P T 

By substituting into eqs. (J4.21)l . ()4.22|1 . ()4.23j) and projecting one obtains 

[□ + V L + (^ 2) -§0i)A£]0 = O 

[□ + m 2 T + V T ]u h T = , [□ + ml + V L ]ukL = 

V L = § 0i [A£ + J L ] + \ (N - 1) e 5 J T , V T = \9, [A^ + Jl] + [\{N -1)03 + 1 6 2 ] J T 

where ^ = <fi 2 and Jl and Jt have the obvious meaning [see eq. ([4.10)1 ] . Now let us rescale 
the coupling, the background field and the vacuum expectation value as prescribed by the 
standard large iV procedure 

Ao^Ao/iV 4>{x) -> \^(j)(x) v 2 ->Nv 2 

By taking the limit and using the symmetry conditions one can see that the correct large 
N equations 



[□ + V T }<P = , [□ + V T ]u kT = 

d 3 p 



V T = |AA^ + |A / 

J\p\<A 



A 1 — A 1 - | Itt(Qs) 



(2vr) 3 



u p (x) u p (y) 



(4.35) 



are recovered, provided the limits of the free parameters satisfy the following relations 



I v 2 9^(0) -m\ -> 

r(2)i-i 



N (0! - 6 A ) ^ i A [1 - | A (4.36) 
NQf^lXo [l-lAoJ^]- 1 

Suitable conditions on the functions Fa then follow from eqs. ([4.33)1 . Notice that in the 
integrals I^t in eqs. ([4.35)1 and (J4.3fjj) we have m\ = 0, as consistent with the limits of 
the parameters. 

To conclude, in the HF approximation tt\ and the 6 a fulfill positivity conditions. 
More precisely they are real and positive when evaluated at the purely spatial value of 
the momentum q s [see eq. (13.23)) ]. 
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We now proceed in defining our modified HF approximation by making some further 
sensible requirements that are not satisfied by the HF approximation. 



First of all we require our approximation to be gapless. That is to say we require that 
the transverse mass of the internal propagator is zero. As a consequence by eq. (|4.28j) we 
have 

ml = \ v 2 4 (O) => F 5 (x, y, z;N) = \ F 4 (x, y, z; N) (4.37) 
Due to this condition the transverse internal propagators are now massless and therefore 

(2) 

the Goldstone bosons loop integral J^j. is logarithmically IR-divergent. The symmetry 
relations in eq. ()4.37|) and eq. ()4.31|) define the longitudinal mass m| and in terms of 

(2) 

the zero momentum values of #4 and £l 4 respectively. To avoid IR divergences we then 
require that 64 and do not depend at all on IjQ. Thus we can write 

F 3 (x,y,z;N)=K 3 (x,y;N) , F^x,y, z; N) = K 4 (x,y; N) 

where Ka{x,V] N) (A = 3,4) are arbitrary functions of two variables only (plus N). The 
same argument apply to the parameters //W and of the manifestly symmetric form 
in eq. (|4.32|) : we require that \ Q\ — ^1 and do not depend on 1^,, which in turn 
implies 



F 2 (x,y,z;N) 



K 4 (x,y;N) K 4 (x,y;N) 



1 + \ K,(x, y; N)(y - x) [i+^k^n]^-^ + ^ y _ 



1+2 K 4 (x,y;N)(y~z) 



F 1 (x, y, z; N) = -\ 1 ^1^3 + Kl {x, y; N) 

IN — 1 — 4 h K 2 (x, y: N) 

where, again, KA{x,y; N) (A = 1,2) are arbitrary functions of two variables (and N). 

Now, as in sec. 13.31 we require the renormalizability with the 1-loop beta function. 
That is, we assume the following RG equation for the bare coupling A 

9Xo - N + 8 \ 2 + 0(A->) 

a!oiA-^4^ Ao + 0(A } 

and ask that the free parameter function do not depend on log A. By the same 
procedure of sec. 13.31 and using 

OI LL _ 01 TT /o/a-Ia __ ui TL 1 



+ 0(A- X ) = + O(A^) = -— 1 + OiA- 1 ) 



<91ogA <91ogA v <91ogA v ' 8tt 2 

we obtain the general forms 

K (x vN)= fM^y)l3 ? = 13 4 

A ' l-l(N + 8)x ' 3 ' ' 

K 2 {x,y;N) = f 2 (a(x,y);N) 
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where 

t \ x ~ y 
a(x,y) = r- : — 

(2) (2) 

Recalling that x = and y = XqIj,£ we see that the Fourier transform of a positive 

definite for purely spatial value of the momentum. Notice also that a vanishes when 
v — > since 1^1 and 1^1 coincide in this limit. 

Notice that now, in contrast to the N = 1 case studied in sec. 13.31 [see eqs. ([3.380 ]. the 
requirement of renormalizability does not fix completely the form of the free parameters. 
This occurs because for N > 1 there exist three distinct finite parts associated with the 

(2) (2) (2) 

logarithmic cutoff divergence of I L [, and 1^. So now we are left with four arbitrary 
functions Ja of one variable (and N) to determine. 

It is convenient to rewrite all the constraints previously found in terms of this new 
parametrization. After some straightforward, albeit rather involved calculations, one 
obtains 

• Perturbative matching: 

h(a; N) = i [{N - 1) - 3 £(0; N) - 9 /J(0; N)\ a + 0(a 2 ) 

f 4 (a; N) = 2/3 + 0(a) , / 3 (a; N) = 3{N \ 8) ~ \ a + 0(a 2 ) 

f 2 (a;N) = -(N+l) + 0(a) 

• Matching in the large- N limit: 

h(a/N; N)^ 0, f 4 (a/N; N) - /(a) , f 2 (a/N; N) + N ^ 1 - ^ f( a ) 
where the arbitrary function / enters only the decoupled longitudinal sector. 

(2) 

• Positivity constraints on Q\ and 9 a : on the real interval < £ < 6/(N + 8) 

o</ 4 (C;AO<^, / a (C;iV) + JV-i<o, /i(C;^)>o 

Actually these conditions are slightly stronger than strictly necessary. In any case 
they also guarantee the positivity of all the 6 a parameters (which holds true in the 
usual HF approximation as well). 

One last requirement is that the N — 1 case should be recovered for all the parameters 
that have meaning also in this case. These are Q\ , Q\ and #4 (that, evaluated at zero 
momentum, determines the longitudinal mass according to eq. ([4.370 ). When N — 1 they 
should have the following form [see eqs. ([3.380 ] 

0! = 12 nS 2) = § 6> 4 = A [1 - f A /fi]- 1 



that is, in terms of the functions /'s, 



/ 4 (a;l) = f , / 2 ( a; l) = -l— L , h(a;l) 



2 -la . , a 



3 1 — a 9 — Sa 
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Notice also that all the possible choices of free parameters coincide when v — > 0, since 
a vanishes in this limit and the perturbative constraints uniquely fix the form of the /'s 
when a = 0. Moreover one can verify that the resulting v = improved HF coincide 
with the massless limit of the improved approximation in the unbroken symmetry phase 
defined in ref. |3D] . 

In conclusion all possible forms of the functions /'s that fulfill the above matching 
constraints define improved Hartree-Fock resummations with the required features of 
gaplessness and renormalizability. One simple choice for the /'s is 

1 - ±a 

/ 4 (a;JV) = 2/3 , / 2 (a; N) = -(N + 1) — + (N - 1) a 

2 1 ( 438 ) 

h{a]N) = 9-M ' /3(a;iV) = 3(ivT^)IT^ 

As already remarked, this choice is not unique. For example we can consider a second 
form 

9/3 1 — ±rv 14- Af + 8 rv 

1 4- — — ^ f> 1 — ol 1 4- ^ rv 

6 6 (4 391 

a „ , _ 2 1 v ' ; 



f 1 (a;N) = — — - , / 3 (a;iV) 



9 -2a ' JJV ' y 3(iV + 8) 1 + 



Af+8 



a 
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As already explained in sec. 13.31 our modified resummation adds leading logarithm con- 
tributions of diagrams that are not present in the usual HF approximation; then the 
two forms just provided of the free parameters, as well as all the other possible ones, 
correspond to different choices of the associated finite parts. 

Given one specific choice for the free parameters, we can proceed in applying the 
coupling constant renormalization condition in eq. (j4.17|) . Notice that the consistence of 

(2) 

this renormalization requires that Q\ is monotonically growing with Ao for any given 
purely spatial momentum. We have omitted to include this requirement in the previous 
general discussion since it would lead in general to rather complicated constraints. It holds 
true for the two simple examples given in eqs. ()4.38|) and eqs. ()4.39|) . as one can explicitly 
check. Moreover it holds true in general (i.e. for any improved HF resummation) for 
scales such that s 2 ^> v 2 , since we have in this limit 

Qf\q s ) = - + 0(m 2 /s 2 ) 

where 1^ stands for anyone of 1^1, I^l an< ^ ^-tt- ^ ne ^normalization condition thus 
defines the bare-to-renormalized relation 

A = e A ° . + 0(m 2 /s 2 ) 

which is the usual 1 loop RG-invariant relation up to 0(m 2 /s 2 ) terms. This shows that 
the direct coupling renormalization condition is approximately scale invariant for high 
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renormalization scales. To obtain complete scale invariance it is enough to slightly modify 
the bare-to-renormalized parameterization by changing the renormalization condition in 
the following way 

6Q?\q s ) = 6\ F 1 (\ lf ) L (q s ), X I^ L (q s ), A /$(?-); N ) 

r ei 6\F 1 (\J%(q s ), A/gfe), \j? T (q s ); N) (4-40) 
= X + 0(m 2 /s 2 ) 

where 

4?(p) = 4 2 2-4 2 2fe)i, 4 2 2(p) = 4 2 2-42fe)i, 4?b) = 4l-4 2 2te)i 

are the properly subtracted finite loops (notice that the subtraction term is always the 

(2) 

purely longitudinal I L [{q s )). Then, using the parameterization of F\ in terms of the 
fi(a; N), one can verify that the bare-to-renormalized relation reads exactly 

A = ^— (4.41) 

i-^Ao42fe) 

Comparing this with eq. ()4.19|) we can see that this has the correct 1— loop-RG improved 
behaviour with a specific choice of the finite parts. Then all the free parameters take the 
following manifestly finite forms 

Qi = AFi(a4l> a4l) a4tj N) , 0i = AF 2 (a4l) a4l; a4t! N) 
= XFi{\J LL ,XJ TL ,XJ TT ]N) , 64 = X F±(XJ LL , XJ TL , XJ TT ] N) 

Moreover, consistently with the renormalization condition in eq. (|4.40|) . we can define the 
running coupling constant by means of the following equation 

X(p)F 1 (X(p)J L 2 l(p), A( P )4 2 2(p), A(p)4 2 )(p); N) = n?(p) 

whose solution is simply the extension of eq. (J4.41)) to arbitrary momentum 

Hp) 



i-^ajS( p ) 

In term of this and of the subtracted integrals we can write the renormalized form of the 
Fourier transform of a: 

a(p) = X(p)[j! r 2 l(p)-J L %)} 

Notice that a becomes small for large (but smaller than the Landau pole) spatial mo- 
mentum q s . More precisely it is of order 0(X(q s ) v 2 /s 2 ). Because of the Landau pole, the 
cut-off of the theory cannot be removed but should be fixed to values suitably smaller 
than the pole (i.e such that X(q\) is of order one). The mass scale of the theory, that 
is v, is much smaller than the cut-off itself. Therefore for spatial momenta with values 
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near the cut-off a is small. Then the perturbative (in a) matching conditions assure 
that the different allowed choices of the free parameters give the same results to order 
0(X(q\) v 2 /A 2 ). In this sense we can say that all the class of improved HF approximations 
shares the same UV behaviour. 

One can easily check that if we renormalize at scale s with coupling constant A and 
at scale s' with coupling constant A' = \(q s >) we define the same bare coupling constant. 
That is to say that the bare-to-renormalized relation in eq. (j4.41j) is RG-invariant. As a 
consequence the expressions of the parameters in terms of Ao and 1^ in eqs. (14.33)1 can 
be thought as manifestly scale invariant definitions. 

5 Conclusions 

In this paper we extended to the case of spontaneously broken symmetry the improvement 
of the HF approach to the O(N) scalar theory that was proposed in ref. [301 - m contrast 
to the standard HF approximation, our improved one is renormalizable, RG-invariant and 
correctly gapless in the Goldstone sector. However, it is not unique, except for N = 1, 
because the mass difference between the longitudinal sector and the transverse Goldstone 
sector allows for a richer structure of renormalization finite parts that cannot be fully 
restricted by the requirements of renormalizability and RG-invariance. As a consequence, 
an entire class of improvements is identified, parametrized by the four functions /a (a; N) 
which must satisfy some further constraints explicitly written in the previous section. 
It is important to stress that, albeit of functional type, the remaining freedom is much 
smaller than what simple dimensional analysis would allow even when simple momentum- 
independent observables are concerned (the momentum dependences in our approximation 
are anyway fixed by construction to be of mean-field type). This makes it difficult, 
although not a priori impossible, to identify consistent schemes to further constrain the 
functions /a(o;; N) by requiring agreement in suitable calculations beyond one loop. Work 
on this direction is in progress. 
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